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Abstract. This paper analyzes the reconstruction of diffusion and absorption 
parameters in an elliptic equation from knowledge of internal data. In the application 
of photoacoustics, the internal data are the amount of thermal energy deposited by high 
frequency radiation propagating inside a domain of interest. These data are obtained 
by solving an inverse wave equation, which is well-studied in the literature. 

We show that knowledge of two internal data based on well-chosen boundary 
conditions uniquely determines two constitutive parameters in diffusion and 
Schrodinger equations. Stability of the reconstruction is guaranteed under additional 
geometric constraints of strict convexity. No geometric constraints are necessary when 
2n internal data for well-chosen boundary conditions are available, where n is spatial 
dimension. The set of well-chosen boundary conditions is characterized in terms of 
appropriate complex geometrical optics (CGO) solutions. 

Keywords. Photoacoustics, Optoacoustics, diffusion equation, inverse problems, 
internal data, stability estimates, complex geometrical optics (CGO) solutions. 

1. Introduction 

Photoacoustic tomography (PAT) is a recent hybrid medical imaging modality that 
combines the high resolution of acoustic waves with the large contrast of optical waves. 
When a body is exposed to short pulse radiation, typically emitted in the near infra-red 
region in PAT, it absorbs energy and expands thermo-elastically by a very small amount; 
this is the photoacoustic effect. Such an expansion is sufficient to emit acoustic pulses, 
which travel back to the boundary of the domain of interest where they are measured. 

A first step in PAT is therefore to reconstruct the amount of deposited energy from 
time-dependent boundary measurement of acoustic signals. Acoustic signals propagate 
in fairly homogeneous domains as the sound speed varies very little from one tissue to the 
next. The reconstruction of the amount of deposited energy is therefore quite accurate 
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in many practical settings. For references on the practical and theoretical aspects of 
PAT, we refer the reader to e.g. P[a[6l[7l[9l[8l[T0l[ia[l3l[Ta[l5l[T6l[I71[2^ 

Once the amount of deposited energy has been reconstructed, a second step consists 
of inferring the optical properties in the body. Although this second step is less studied 
mathematically, it has received significant attention in the biomedical literature; see e.g. 
[HllHllT]. The reconstruction of the optical parameters is very useful in practice because 
the attenuation properties of healthy and unhealthy tissues are extremely different 
[3 [20]. The combination of high spatial resolution of the acoustic signals and large 
contrast of the optical parameters allows us e.g. to image blood vessels at high resolution, 
which is in turn important for assessing the status of cancerous tissues. 

Near infra-red light is best modeled by radiative transfer equations. What we can 
reconstruct on the optical coefficients, seen as constitutive parameters in the radiative 
transfer equation, from measurements of heat deposition for all possible illuminations 
of the domain has recently been analyzed in [1]. The latter paper considers the 
continuous-illumination setting. In other words, heat deposition is measured for all 
possible incoming radiation condition, which can be controlled at the boundary of the 
domain. 

In this paper, we consider the diffusion approximation to radiative transfer, which 
is accurate for radiation propagation in highly scattering media. Such an approximation 
is typically valid for propagations of radiation over one centimeter or more [2], [3] . In such 
a simplified setting, the unknown optical parameters are the spatially varying diffusion 
and attenuation coefficients. Provided that the diffusion coefficient is known at the 
domain's boundary, we show that these two coefficients are uniquely determined by two 
well-chosen illuminations at the boundary of the domain. An explicit reconstruction 
procedure is presented, which involves solving a first-order equation with a vector field 
explicitly constructed from the internal data (the amount of deposited energy). 

The stability of the reconstruction of the optical coefficients from two internal data 
is established under geometric conditions of strict convexity on the domain of interest. 
In the presence of 2n well-chosen boundary conditions, where n is spatial dimension, 
such geometric constraints can be removed and constructions are possible on essentially 
arbitrary domains of interest with well-defined boundary. 

Mathematically, the inverse problem is an inverse diffusion problem with two 
internal measurements. By using the standard Liouville change of variables, the diffusion 
equation is replaced by a Schrodinger equation with unknown potential and with internal 
measurements involving a second unknown function. By adapting the theory of complex 
geometrical optics solutions to this setting, we are able to obtain uniqueness and stability 
results for the inverse Schrodinger problem. By means of the inverse Liouville change of 
variables, we are able to conclude on the uniqueness and stability of the reconstruction 
of the optical parameters for the photoacoustic tomography problem in the diffusive 
regime. 

The rest of the paper is structured as follows. Section [2] presents the photoacoustic 
tomography problem and the main results we obtain in this paper. The inverse 
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Schrodinger problem and the explicit reconstruction algorithms are addressed in section 
[3l How the latter results are used to solve the inverse diffusion problem is explained in 
section HI 



2. Photoacoustic tomography and main results 

The propagation of radiation in scattering media is modeled by the following diffusion 
equation 

1 d 

-—u - V ■ Dix)Vu + aa{x)u = 0, a; G X C M", t > 

c at [1) 

u = g X e dX, t > 

where X is an open, bounded, connected domain in with boundary dX (embedded 
in M"), where n spatial dimension; c is light speed in tissues; D{x) is a (scalar) diffusion 
coefficient; and cra(x) is an attenuation coefficient. The source of incoming radiation 
is prescribed by g{t,x) on the boundary dX and is assumed to be a very short pulse 
supported on an interval of time (0,?]) with crj of order 0(1). The amount of energy 
deposited is proportional to attenuation and is given by 

H{t,x) = aa{x)u(t,x). 

A thermal expansion (assumed to be proportional to H) results and emits acoustic 
waves. Such waves are modeled by 

with Cg the sound speed and /5 a coupling coefficient assumed to be constant and known. 
The pressure (potential) p{t, x) is then measured on dX as a function of time. Using 
as in [1] the difference of time scales Cg <^ c, which shows that radiation propagation 
occurs at a much faster time scale than acoustic wave propagation, we can show that 

H{t,x) ^ Ho{x)6o{t), Ho{x) = aa{x) / u{t,x)dt. 

Jr+ 

We now have a well-posed inverse wave problem, where measurement of n— dimensional 
information p{t, x) for t > and x G dX allows us to reconstruct the n— dimensional 
spatial map Hq{x). We assume this step, which has been analyzed in great detail in the 
literature, completed; see e.g. P fTH fT5l fTTl [21] . 

Denoting by g{x) = J^^ g(t, x)dt and u{x) = J^^ u(t, x)dt, we thus observe that the 
photoacoustic problem in the diffusive regime amounts to reconstructing {D{x),aa{x)) 
from knowledge of 

Hj{x) = aa{x)uj{x), 1 < j < J, (3) 
for J G N* illumination maps gj{x), where Uj is the solution to the steady-state equation 

-V ■ D{x)Vuj + aa{x)uj = 0, x G X C M'', 
Uj = gj X G dX. 
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The set of internal data is given by: 

d = {dj)i<j<j, dj{x) = aa{x)uj{x). (5) 

The main problem of interest in this paper is the uniqueness and the stability of 
the reconstruction for the Inverse Diffusion problem with Internal Data: 

ISID: Reconstruction of {D{x),aa{x)) from knowledge of d = ('ij)i<j<j on X for a 
fixed collection of illuminations g = {gj)i<j<j prescribed on dX. 

We start with the case of two real- valued measurements J = 2. The reconstructions are 
based on the construction of vector fields that are well defined only when the optical 
coefficients are sufficiently smooth. More precisely, let k > 1 and assume that 

/D G r = i/t+fc+2+^(x) c C^'+^X), aa e C''+\X), e>0. 

Let then t G M" be a constant vector with \t\ sufficiently large. 

Let g = ((71, (72) be a given illumination and d = di + id2 = aa{x){ui + iu2) the 
corresponding internal data, where the real-valued solutions Uj solve (jl]) with boundary 
conditions gj for j = 1,2. We introduce the vector field and scalar quantity 

f3 = ^ ^ _-^/\^y (6) 

Here, d = di — id2 is the complex conjugate of 0. Then we will see in section |3] that for 
jjL = D~^aa, we have 

/?-V/i + 7/i = 0. (7) 

Since /i is known on dX as D and cTq = dj/gj are known on dX, the above equation is a 
well-posed equation for fi. When the integral curves of the vector field /3 map any point 
inside a; to a point xo{x) G dX, then ([7j) uniquely characterizes /i inside the domain. 

Let V = ^/D^u and define q = —v~^Av, which we will show is well-defined. Then 
it turns out that 

-A^/D - qVD = fi, xeX. 

This provides an elliptic equation for \/D that admits a unique solution since y/D is 
known on dX. 

The principal difficulty is to ensure that (3 is constructed in such a way that ([71) 
allows for a unique solution /i. We will show that there exists an open set of boundary 
conditions g on dX that allows us to do so. We define the set of coefficients {D, aa) G Ai 
as: 

M = {iD,aa) such that (/D,crJ eY x C''^\X), \\^/D\\y + \\aa\\c>=+^x) < M}. (8) 

The main results for the Inverse Diffusion problem with Internal Data (IDID) are 
then as follows. We start with a uniqueness results: 
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Theorem 2.1 Let X he an open, hounded, domain with C"^ houndary dX . Assume that 
{D,aa) and {D,da) are in M. with D\gx = D\dx- Let d and d he the internal data in 
([5]) for the coefficients {D,aa) and {D,aa), respectively and with houndary conditions 
((?j)j=i,2- Then there is an open set of illuminations g G {C^'°'{dX)y for some « > | 
such that if d = d, then {D,aa) = {D,aa)- 

The above resuh shows uniqueness of the reconstruction of the optical coefficients 
but does not imply stability. Without additional geometric information about dX, the 
above procedure may yield unstable reconstructions. However, provided that 2n well- 
chosen measurements dj for 1 < j < J = 2n are available, then n vector fields similar 
to (3 above may be constructed in order to form locally a basis of vectors in M". In such 
a setting, the following stability result holds. 

Theorem 2.2 Let k > 2 and let X he an arhitrary hounded domain with houndary 
dX of class C^'^^. Assume that {D,aa) and {D,aa) are in M. with D\qx = L)\qx- 
Let d = {di, . . . ,d2n) and d = {di, . . . ,d2n) he the internal data constructed in ([5]) 
for the coefficients {D,aa) and {D,aa), respectively and with houndary conditions 
g = {gj)i<j<2n- Then there is an open set of illuminations g G (C''^'°(9X))^" and a 
constant C such that 

\\D - D\\ck(^x) + Ika - (^a\\c''{X) < C\\d - (cfc+i(x))2" • (9) 

Stability is therefore ensured without geometric constraints provided that enough 
measurements are available. When only two measurements are available and are of the 
form described in Theorem 12.11 we can still get a stability result under the following 
geometric hypothesis: 

Hypothesis 2.3 There exists R < oo such that for each xq G dX , which we assume is 
of class C"^ , we have X C Bxq{R) where Bxq{R) is a hall of radius R that is tangent to 
dX at xo G dX. 

Then we can show the following result. 

Theorem 2.4 Let k > 3. Let X satisfy Hypothesis \2.3\ with houndary dX of class C^^^. 
Assume that {D,aa) and {D,aa) are in M. with D\qx = D\dx- Let d and d he internal 
data as ahove for the coefficients {D,aa) and {D,aa), respectively and with houndary 
conditions g = {gj)j=i,2- Then there is an open set of illuminations g G (C'^'"((9X))^ 
and a constant C such that 

\\D - D\\ck~i(^x) + \Wa - O-aWck-^X) < C\\d - rf||(cfc(X))2- (10) 

The above three theorems are proved in the following two sections. They show 
that the internal data d for illuminations g with J > 2 provide stable reconstructions 
of the optical coefficients under the assumption that the illuminations are well-chosen. 
How these illuminations are chosen will become more explicit in the next two sections. 
The characterization of the open set of illuminations is however not very precise. The 
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main features of the result are as follows. For coefficients in M.^ there is a minimal 
value of |6| that ensures that we can construct a vector field (3 with Property P, which 
means its integral curves map any point in X to an point in dX. Such a vector field 
is constructed by means of complex geometric optics solutions, which depend on the 
unknown optical parameters. The illuminations must then be chosen sufficiently close 
to the trace on dX of the above vector field to ensure that they generate another vector 
field with Property P. Closedness is therefore not characterized in very explicit means. 
It remains an interesting question to obtain a priori constraints on the illumination that 
will ensure that the resulting vector field satisfies Property P. 

In the next section, we consider a similar problem for the Schrodinger equation. 
How the latter results are used to prove the theorems stated above is described in 
section HI 



3. Inverse Schrodinger with Internal Data 

Let X be an open, bounded, connected, domain in M*^, where n is spatial dimension, 
with smooth boundary dX. We consider the Schrodinger equations 

Am,- + quj = Q X 

(11) 

Uj = Qj dX, 

for 1 < J < J. Here, J G N* is the number of illuminations and q is an unknown 
potential. We assume that the homogeneous problem with gj = admits the unique 
solution M = so that A = is not in the spectrum of A + g. We assume that q on 
X is the restriction to X of a function q compactly supported on and such that 
q G ift"'''^"'"^(]R") with e > for /c > 1. Moreover we assume that the extension is chosen 
so that 

WqWh^+^+^x) ^ C{X,k,n)\\q\\j^^+k+e^^„y (12) 

for some constant C{X,k,n) independent of q. That such a constant exists may be 
found e.g. in fiSl Chapter VI, Theorem 5]. 

We assume that gj G C^''°(X) with a > \ and dX is of class C^'^^ so that (ITTIl 
admits a unique solution Uj G C^^^{X) p!Tl Theorem 6.19]. The internal data are of 
the form 

dj{x) = fi{x)uj{x), X, I < j < J- (13) 

Here /i G C'^+^(X) verifies < yUo < /i(x) < /ig ^ for a.e. x E X. 

The inverse Schrodinger problem with internal data (ISID) consists of reconstruct- 
ing (g, /i) in X from knowledge ofd= {di, . . . , dj) G (^C^'^^{X)Y for a given illumination 
g = igj)i<j<j- We will mostly be concerned with the case J = 2 and J = 2n with gj, 
and hence dj real-valued measurements. 
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3.1. Complex Geometrical Optics Solutions 

The analysis of ISID carried out in this paper is based on the construction of complex 
geometrical optics (CGO) solutions. When g = 0, CGOs are harmonic solutions of the 
form e^'^ for p G C" such that p- p = 0. When q ^ 0, CGOs are solutions of the following 
problem 

Aup + qup = 0, Up e^'^ as |a;| oo. (14) 

More precisely, we say that Up is a solution of the above equation with p ■ p = and the 
proper behavior at infinity when it is written as 



u 



,{x) = eP--{l + ijp{x)), (15) 



for ipp & L'j a. weak solution of 

A^p + 2p-V4!p = -qil + ^p). (16) 

The space for 5 G M is defined as the completion of C^(]R") with respect to the norm 
II ■ \\i2 defined as 



\25\ 



u\\q={l {xy'\u\'dx), {x) = {1 + \x\')^ . (17) 



Let — 1 < (5 < and q G L'j^-^ and {x)q G L°°. One of the main results in [IS] is that 
there exists rj = ri{6) such that the above problem admits a unique solution with ip E L'j 
provided that 

\\{x)q\\Loo + 1 < ?7|p|. 

Moreover, HV^Uli < C|p|^^ ||g||^2^^ for some C = C{6). In the analysis of ISID, we need 
smoother CGOs than what was recalled above. We introduce the spaces if| for s > 
as the completion of C^(M") with respect to the norm || ■ defined as 



11^11^^!= / {xy'\{I-A)-2u\'dx]\ (18) 

Here (/ — A)5m is defined as the inverse Fourier transform of (O^'^iOy where u{^) is the 
Fourier transform of u{x). Then we have the following 

Proposition 3.1 Let —1<6<0 and k E W. Let q G " and hence in H^^^ 

and p be such that 

||g||^t+.+. + i<^|p|. (19) 



Then tpp the unique solution to (fT6l) belongs to "'"'^"''^ and 



\p\\\iJp\\^^+,+.<C\\q\\^^,+., (20) 



for a constant C that depends on 6 and rj. 



Inverse Diffusion Theory of Photoacoustics 8 

Proof We recall [IH] that for |p| > c > and / G L^^-^ with -1 < 5 < 0, the 
equation 

(A + 2p-V)7/; = / (21) 
admits a unique weak solution tp & Lj with 

ii^iu. <c(5,c)|pr^ii/iu.^^. 

Now since (A+2p- V) and (/—A)* are constant coefficient operators and hence commute, 
we deduce that when / G Hs+i ^'^y > 0, then 

||^||h|<C(5,c)|p|-^||/||h|,,. (22) 
The solution to (fT6l) is known to admit the decomposition 

oo 

with = 1. Let s = ^ + k + e. Assume qipj-i G -f^l+i which is true for j = by 
assumption on q. Then ipj G Since H'^ is an algebra, we want to prove that 

\\qi^j\\Hi^^ < \\q\\H^Uj\\Hi- (23) 

Indeed, decompose into cubes. On each cube B, (x)^* is more or less constant up to 
a C^^*. Now H^{B) is an algebra so that < Since (x) is 

more or less constant and equal to (x^), 

(a:B)^'+^||g«||^.(B)<C||(x)(/-A)fg||i.(^)||(x)^(/-A)i«||i.(^) 

It remains to sum over all the cubes B to get the result. When the size of the cubes 
tends to 0, the constant C tends to 1, which yields (l23l) . This and (l22l) show that 

< ciprigii^f+i+.ii^/'i-iiii/|- 

By selecting rj such that C|p|~^||g|| ^+k+E < |, we obtain 

It remains to sum the geometric series to obtain the result. □ 

We now want to obtain estimates for %l)p and Up restricted to X. We have the following 
result. 

Corollary 3.2 Let us assume the regularity hypotheses of the previous proposition. 
Then we find that 
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Proof. On the bounded domain X, (x) is bounded above and below by positive 
constants. Since q is compactly supported on M", we obtain thanks to (fT2|) that 

|p|||^plli^5+fe+-(x) - C'll'?ll//5+'--+^(R") - 

Now we have 

AV^p = -2p- VV'p-g(l+V'p)• 
By elliptic regularity, with X' a smooth domain in such that X C X', we find for 
all s = I + /c + e, that 

< C'l|2p ■ V^p - g(l +^p)||//=-i(x') + II^pIIh^cx')- 

The latter is bounded by |p| ||^/'p||_ff»(x') + Iklln^-M^') ll^pll-'^^^H-^') since s — 1 > | so that 
H^~^{X') is a Banach algebra. By using the above bound on ||^/'p||/i-s(j5s:'), with C{X) 
replaced by the larger C(X'), we get the result. □ 
By Sobolev embedding, we have just proved the: 

Proposition 3.3 Under the hypotheses of Corollary \3.2\. the restriction to X of the 
CGO solution verifies that 

IpIII^pIIch^) + il^pllcfc+i(^) - CM\H^+''+^xy 

(25) 

We recall that q satisfies the constraint (fT9l) . 

We are now in a position to prove the main result of this section. 

Theorem 3.4 Let Up^ for j = 1,2 be CGO solutions with q such that 0191) holds for 
both pj and k > 1 and with Cq'^\pi\ < \p2\ < Co|pi| for some Cq > 0. Then we have 

/3 := J-e-'-^^^^-y^UyUp, - Up^Vup,) = ^^-^ + K (26) 

Z\Pl\ \ / A\pi\ 

where the vector field h satisfies the constraint 

WHc^ix) < ^, (27) 

IpiI 

for some constant Cq independent of pi^2- 
Proof. Some algebra shows that 

h = ^^(^. + i-^ + i-M + ^''"'^ * " r^'^ (28) 

^iPlI ^iPlI 

We know from Proposition 13.31 that |pj||'?ApJ and jV-i/^pJ are bounded in C^{X) for 
j = 1, 2. This concludes the proof of the theorem. □ 



Inverse Diffusion Theory of Photoacoustics 



10 



3.2. Construction of vector fields and uniqueness result 

Let us consider two internal complex- valued data di^2{x) obtained as follows. We assume 
that we can impose the complex-valued boundary conditions gi^2 £ C^'"((9X;C) and 
define the solution mi 2 of 

Auj + quj = 0, X, Uj = gj, dX, j = 1,2. (29) 



Note that the real and imaginary parts of ui^2 may be solved independently since ( 1291) is 
a linear equation. We then assume that we have access to the complex-valued internal 
data dj = fiUj on X for j = 1, 2, where Ui^2 are the solutions of fl29|) with boundary 
conditions gi^2- We recall that fj, G and is bounded above an below by positive 

constants. We verify that 

uiAu2 — U2AU1 = 0. 

Introducing u = -, which is well defined since /i is bounded away from 0, and using 
Uj = udj, we obtain that 

2{diVd2 - d2Vdi) ■ Vv + {diM2 - d2Mi)v = 0. 

This is equivalent to 

/3d ■ V/i + 7d/i = 0, (30) 

where 

/3d := xix){diWd2- d2Vdi) 

7d := 7;X{x>{d2Adi - diArfa) = . 

Here, x{^) is a smooth known complex-valued function with |x(3;)| uniformly bounded 
from below by a positive constant on X. Note that by assumption on /i, we have that 
/3d G {C\X- C)T and 7d G C\X- C). 

A methodology for the reconstruction of (/i, q) is therefore as follows: we first 
reconstruct /i using the real part or the imaginary part of (l30l) for then 3?/^d and 53/^d 
are real- valued vector fields since fi = d/g is known on dX. When /x is reconstructed, this 
gives us explicit reconstructions for M12 = 1^1,2//^ and we may then reconstruct q from 
the Schrodinger equation. Such a method provides a unique reconstruction provided 
that the integral curves of (the real part or the imaginary part of) f3d join any point in 
X to a point Xq{x) G dX, where fi is known. We thus need the vector field $ to satisfy 
such properties. CGO solutions will allow us to construct families of vector fields $d 
with the required properties. 

Let us consider two CGOs jWith parameters pi^2- Let di 2 be the complex-valued 
corresponding internal data. Let us decompose as before 

u.^ix) = e'^r^il + ij.^ix)), V«p,(x) = e^^-^((l + z^pjp, + Vz^p,). 

Let us choose x(x) = e'^^^'^^'^^'^ in (13T]) . Then we find after some algebra that Pd in (l3T!l 
is given by 

/3d = /i'((pi - P2)(l + V^pJ(l + V^pJ + Vz/.p,(l + zApJ - Vz/'p,(l + z/^pj). (32) 
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We may then define 

/3:=-^/3, = /i2^^ + A 7:= 7^7., (33) 
2|pi| 2|pi| 2|pi| 

where h is defined as in (l28l) . Then, we deduce from Theorem [33] that is bounded 

uniformly in C^{X\ C). When |pi| is sufficiently large, then /? is close to /^^ %\p[\ ' '^hich 
is a non- vanishing vector when pi 7^ p2. Provided that the real part or the imaginary 
part of $ does not vanish, then (l30ll gives an equation for p that can be uniquely solved 
since p is known on dX. 

Note that the data are complex valued. The only possibility to construct two 
different complex valued data with two real valued data is to assume that 62 = di, the 
complex conjugate of di. For the construction of CGOs, this implies that we choose 
p2 = 7h- Indeed, we verify that W^ = Up since for p = t + it-^ with \t\ = and i-i^ = 0, 
we have eP'^ = e^'^ and ipp = tpp by uniqueness of the solution to the equation satisfied 
by ipp. This implies then that /3 defined in f l55]l with p2 = pi is given by 

pp = ifxH^ + p2/i. (34) 

As soon as \p\ > Co so that ||/i||c;o(X) < 1, we obtain that any point in X is connected 
to a point in dX by an integral curve of /3p := Q$p. 

Note that Up solves (129|) with the unknown boundary condition Up\Qx G C'^'^^dX; C) 
for some a > | since Up is known to be a little more regular than being of class 
C^'+^(X; C) by construction (since e > 0). 

Let us now define boundary conditions g G C'^'"((9X; C) such that 

\\9 - Up\dx\\c''''^idX;C) < ^5 (35) 

for some e > sufficiently small. Let u be the solution of ( l29l) with g as in ( l35l) . By 
elliptic regularity, we thus have 

\\u-Up\\ck+nx;C) <Ce, (36) 

for some positive constant C. Define the complex valued internal data d = fiu. Since 
p G C''+^(X), we deduce that 

IM ~ dp\\ck+i(^x.c) ^ C'oe, (37) 

for Co > 0. Once (i is constructed, define di = d and d2 = d and define /3d and p^ as 
in flHTl) with x(3;) = e~^*'^ and the normalized quantities $ and 7 as in fl33p . Note that 
x(a:) is positive and bounded on X. 
Let us define 

:= = 7 := = ^5>7d. (38) 
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Thanks to (13 7p and (1341) . we obtain the error estimate 

||/3-/i'e||cHx)<C^. (39) 

As a consequence, as soon as \V\ is sufficiently large and e sufficiently small, we obtain 
that /3 ■ > C > so that any point x G X is mapped to a point in dX in a time less 
than |(^|~^diam(X) by an integral curve of /?. 

Moreover, we have the equation with real-valued coefficients: 

/? ■ V/i + 7/i = 0. (40) 

Since ^ = d/g is known on dX, this equation provides a unique reconstruction for /i. 
Let us define the set of parameters 

V = \{fi,q) e C'+^X) X iJt+'=+=(X); not an eigenvalue of A + g, 

„ „ „ „ ^ (41) 

\\h\c^+\x) + \\(1\\h^+>^+-(x) < ooj- 

The above construction of the vector field allows us to obtain the following uniqueness 
result. 

Theorem 3.5 Let X be a bounded, open subset o/M" with boundary dX of class . 
Let (/i, q) and {ft, q) be two elements in V . Let 6 G with \t\ > |fio| o,n'd |6o| sufficiently 
large and define p = i + H^ so that p- p = 0. Let Up be the corresponding CGO for q and 
u constructed as above with d = fiu and with e sufficiently small. Let d be constructed 
similarly with the parameters {p,q). 

Then d = d implies that {p, q) = {p, q) . 

Proof. Since the two measurements d = d, we have that p and p solve the same 
equation fj40|) . Since p = p = d/g on dX, we deduce that p = p since the integral 
curves of f3 map any point x & X to the boundary dX. More precisely, consider the 
flow (px(t) associated to j3, i.e., the solution to 

0x(t) = Pi^xit)), = X G X. (42) 

By the Picard-Lindelof theorem, the above equations admit unique solutions since (3 is 
of class C^. And by hypothesis on j3 since (3 ■ i > ( > for |6| sufficiently large, any 
point X is mapped to two points (for positive and negative values of t) on dX by the 
flow (fx{t) in a time less than C~^diam(X). For x G X, let us define x±{x) G dX and 
±t±{x) > such that 

^x{t±{x)) = x±{x) G dX. (43) 
Then by the method of characteristics, p{x) solution of (l40l) is given by 



p{x) = /io(x±(x))e-^o*''''^(^-(^»'^^ 



(44) 
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The solution jl{x) is given by the same formula since V9x(^) = '^x{^) so that fi = fi. This 
implies that u = u since d = d. It remains to use the equation for u to deduce that 
g = g on the domain where m 7^ 0. By unique continuation, u cannot vanish on an open 
set in X different from the empty set for otherwise u vanishes everywhere and this is 
impossible to satisfy the boundary conditions. This shows that the set F C X where 
\u\ > is open and F = X since the complement of F has to be empty. By continuity, 
this shows that q is known on X. □ 

The above result shows that there exists an open set of boundary conditions g close 
to Up\Qx so that data di = d and d2 = d obtained from one complex-valued solution 
u or equivalently from two real valued solutions and '^u, uniquely determine the 
parameters {fi,q)- A more explicit characterization of the open set of illuminations is 
lacking. However, we observe that larger values of q require larger values of |6| in order 
to straighten the vector field /?. Although (l35l) seems to be independent of 1^ and p, in 
fact Up itself grows exponentially with |6| = |p| so that g has to be in the e vicinity 
of an exponentially growing quantity. This means that \t\ has to be sufficiently large 
that the field (3 is sufficiently fiat while at the same time not so large that the imposed 
illuminations become physically infeasible. 

The above uniqueness result does not guaranty stability in the reconstruction. We 
easily verify that the construction provides stability of the reconstruction of fi in most 
of the domain X. However, small changes in the data may generate small changes in the 
field p. This in turn may significantly modify the value of the reconstructed function /i 
at points where (3 is "almost" tangent to the boundary dX. We will see below that under 
some geometric constraints of sufficient convexity of X, the above procedure provides 
a stable reconstruction of the parameters (/x, g). When such conditions are not met, 
we can still obtain stability by acquiring more measurements. Indeed, if a sufficient 
number of vector fields /5 can be constructed at every point so that the span of these 
vector fields is exactly M", then we face a significantly more favorable situation. We 
now consider such a case where 2n real-valued measurements are available. Later, we 
will derive stability results in the two-measurement setting under additional geometric 
constraints. 

3.3. ISID with 2n real-valued internal data 

Let us consider first the setting in which we can access 2n real-valued internal data 
viewed as n complex-valued internal data (since the measurements are linear in u, we 
can measure the real and imaginary parts separately). 

Let us define = \i\ej where (ei, . . . , e„) is an orthonormal basis. We define the 
complex vectors 



Let Up. be the corresponding CGOs. We choose boundary conditions gj such that 



p. = ^1 + iij, 2<j^n, 



Pi 



-ti - it2 = -P2- 



(45) 




(46) 
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for e sufficiently small. We define Uj as the solutions to (12^ with boundary conditions 
Qj. These are n complex- valued solutions whose real and imaginary parts consist of 2n 
real- valued solutions. For 1 < j < n, we define dj = ^uj. We now construct the n vector 
field f3j. For 2 < j < n, the real- valued vector fields and scalar terms are constructed 
as in the preceding section; for j = 1, the vector field is constructed by using p2 = —pi'- 

Pi = ^3ft(d2Vrfi-rfiVd2), 71 = ^3fJ(ciiArf2-rf2Arfi), 

g-2ti-x _ _ ^-2h-x _ _ (47) 

Pj = -^^'^{dj^dj-dj^dj)^ 7j = -j^'^{dj^dj-djAdj), 
for 2 < j < n. As in the preceding section, we verify that 

m-fi%\\cHx)<C^. (48) 

For 1^1 sufficiently large, and thanks to the bound /ig ^ < /i < Po^ we obtain that 
at each point x G X, the vectors Pj{x) form a basis. Moreover, the matrix a^j such 
that Pj = 'Y^ajkCk is an invertible matrix with inverse of class C^{X). In other words, 
we have constructed a vector- valued function r(x) G {C^{X)Y such that ( l40i) may be 
recast as 

V/i + r(a;)/i = 0. (49) 

Finally, the construction of F is stable under small perturbations in the data dj. Indeed, 
invertibility of aj^ is ensured for vector fields close to Pj. Let F and F be two vector 
fields constructed from knowledge of two sets of internal data d = {dj, 1 < j < n} and 
d = {dj, 1 < j < n}. Then we find that 

l|r - r||(cfc(x))n < C\\d - rf||(c"=+i(J?;C))"; (50) 

provided the right-hand side is sufficiently small. 

Let us now assume that X is connected (otherwise, the method applies to each 
connected component) and /i is known and equal to po = d/g for some point xq G dX. 
In other words, we want to solve the over-determined problem 

V/i + F(x)/i = 0, X, /i(xo) = fJ^oixo), Xq G dX. (51) 

Let X G X be an arbitrary point and assume that X is bounded and connected and dX 
is smooth. Then we find a smooth curve that links x to the point xq G dX. Restricted 
to this curve, (!5T!) becomes a stable ordinary differential equation. The solution of the 
ordinary differential equation is then stable with respect to modifications in F (the curve 
between x and Xq is kept constant). The solution p then clearly inherits the smoothness 
of F(x) directly from (15T]) . Moreover, since po{xo) — /io(xo) (with /iq = d/g on dX) 
is small and equation ( 150|) is stable with respect to changes in the value of yUo(xo), we 
deduce that the reconstruction of /i is stable with respect to perturbations in d. 
We may thus state the main result of this section: 
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ll/i — fi 



Theorem 3.6 Let k > 1. We assume that we have access to n well-chosen complex- 
valued measurements and that (/i, q) and {fi, q) are elements in V . Under the hypotheses 
outlined above, and provided that \\up.\dx — gj\\c^,<^{X;C) ^■5 sufficiently small, then we have 
the following stability result: 

'llc"=(x) + Ik - ^llc'=-2{x) < C\\d - ci||(cfe+i(^))2n, (52) 

Proof The inequahty for /x — /i is a direct consequence of the resuhs proved above. 
This provides a stabihty resuh for u = fi~^ and for Uj = udj from the data dj. We thus 
have a stabihty resuh for Auj = —ujq and hence the above stabihty resuh for Uj{q — q) 
since {uj — Uj)q is smaU. 

Now, Up = e'''^{l-\-ipp) does not vanish on X when \p\ is sufficiently large since \p\ipp 
is bounded. When the boundary condition gj — Up.\Qx is small, then by the maximum 
principle, uj does not vanish on X either. This means that either its real part or its 
imaginary part does not vanish everywhere in X. This provides control of g — g in X 
as given in (1521) . □ 

3.4- Vector fields and stability of solutions 

The above construction allows one to stably reconstruct the two functions /i and q 
provided that we have constructed J = 2n well-chosen real-valued boundary conditions 
and collected 2n corresponding internal data. We now return to the reconstruction of /i 
and q in the presence of J = 2 well-chosen real- valued internal data. Such internal data 
are obtained as Theorem 13. 5[ 

We recall that 6 is fixed and p = 6 + it-^. We define u as the solution to (l29l) with 
g close to Up\Qx- The complex- valued internal data are then d = pu. The vector field /3 
and the scalar 7 are then given by 



f] = -^'^(^dVd-dVd), 7 = -——Q{dAd-dAd), (53) 



2|t| ^ ^' ' 4161 

and we verify that 



As earlier, we verify that 



/? ■ V/i + 7/i = X. (54) 



fi't^WcHx) < (55) 



As a consequence, the integral curves of /3 given by ifx{t) map any point x G X to two 
points on dX when |t| is sufficiently large as was mentioned earlier. 

However, the stability of equation (15^ with respect to changes in (3 and 7 is not as 
good as in the presence of n complex internal data. The stability of the reconstruction 
degrades for points x close to xq G dX where Pp{xo) ■n{xo) is close to 0. The stability of 
the reconstruction of fi will however be good when X is a convex domain with "sufficient" 
convexity as established in Hypothesis 12. 3[ We prove the following result: 



Inverse Diffusion Theory of Photoacoustics 16 

Proposition 3.7 Let k > 1. Let fi and fl be solutions of (1541) corresponding to 
coefficients (/3,7) and {/3,'j), respectively, where 

P = fiH^ + ^/i, /3 = fi'l^ + ^~h, 

for h, 7, h, and 7 bounded in C'^(X). 

Let us assume that h\qx = /^o CLiT'd ft\dx = P^o on dX for some functions /io,/io ^ 
C^{dX). Let us assume that X is sufficiently convex so that Hypothesis \2.3\ holds for 
some R < 00. We also assume that \t\ > 60 is sufficiently large. Then there is a constant 
C such that 

ll/i -/illcfe-i(^) < C\\fio\\ck(^dx){\\l3- P\\c^'Hx) + h-l\\c^-Hx)) 

(56) 

+ C'll^O — /^o||c"=(9^)- 

Proof of Proposition 3.7 We recall that ifxit) is the flow defined in (H2|) and that 
x±{x) and t±(x) are defined in ( l43l) . By the method of characteristics, solution of 
f l5^ is given by 

fi{x) = /io(a:±(x))e-^o**'''^(^-(''»'^^ . (57) 

The solution fi{x) is given similarly. We first assume that fiQ = ixq. 
From the equality 

^.{t)-^x{t)= [ [(3{<f,{s))-P{if,{s))]ds, 
Jo 

and using the Lipschitz continuity of (3 and Gronwall's lemma, we thus deduce the 
existence of a constant C such that 

- Mt)\ < ct\\p - pWcoix) 

uniformly in t knowing that all characteristics exit X in finite time and provided that 
(Px{t) and ifx{t) are in X. 

Such estimates are stable with respect to modifications in the initial conditions. Let 
us define W{t) = D^ipxit). Then classically, W solves the equation W = DxP{(px)W 
with W{0) = I and by using Gronwall's lemma once more, we deduce that 



\W - W\{t) < Ct\\DxP - Dx, 



IC^iX), 



for all times provided that (px(t) and (px(t) are in X. As a consequence, since P and /3 
are of class C^{X), then we obtain similarly that: 

and this again for all times provided that (pxit) and (fx{t) are in X. 

However, this does not imply that x+(x) is close to x+(x). When dX is fiat for 
instance, we may very well have that x+{x) is such that n(x+(x)) ■ ipx{t+{x)) is very 
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small and that — is arbitrarily large if /3 is parallel to the surface dX for 

instance. This behavior, however, cannot occur when both (3 and (3 are sufficiently flat, 
which is the case when \t\ is sufficiently large, and when dX is sufficiently curved, which 
is obtained from the existence of i? < oo in Hypothesis 12. 3[ In such a setting, we can 
obtain the following result: 

Lemma 3.8 Let k > 1 and assume that j3 and (3 are C^{X) vector fields that are 
sufficiently flat, i.e., that \t\ is sufficiently large. Let us assume that dX is sufficiently 
convex so that Hypothesis \ 2. 3\ holds for some R < oo. Then we have that 

- + \\t+ - i+Wc^-Hx) < C\\(3 - /3||c7fc-i(x), (58) 

where C is a constant that depends on \i\ and R. 

The above lemma is mostly a consequence of the following result: 
Lemma 3.9 Let Co be the constant defined such that 

|^.(t)| = |V/3(^.(t))/5(^.(t))|<^. 

Let tM be the maximal time spent by any trajectory in X , which we know is bounded. 
Assume that \i\ is sufficiently large that for all x ^ X, 

CqR 1 ^ 

|!| 4|!|2 J\(3{x+{xW ' 

Then we have that for all x & X , 

2R 

'^^^^ - |/?(x^(a:))P(l-p) "^^-^^^^^ ■ ^^^^ 

In other words, the vector field /5(x+(x)) is close to being tangent to dX only when the 
time spent in X is small. 

Let now Xq E X and x G dX and define vq = /3{xq). Assume moreover that 

\x-xo\< CiS"^, \n{x) ■ Vol < C25. (60) 

Then we have 

t+(xo) < C,6, (61) 

for some constant C3 > independent of xq . In other words, a trajectory close to dX 
and almost tangent to dX exits X in a short time. 

We postpone the proof of the above two lemmas to the end of the section. Let us 
conclude the proof of the proposition. We recall that 
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with a similar expression for jl. Since x — > is smooth, by the Leibniz rule it is 
sufficient to prove the stability result for /io(x+(x)) and for /J^''^'' ''y{'^x{.s))ds. It is clear 
from the above lemmas that 

< C'||/^o||c'^(9X)||/5 - f3\\c'^'^(X)- 

Let us now assume without loss of generality that > t+(a;). Then we have 

/ {li'Pxis)) -^{(fxis)))ds = / (7(^^.(s)) -7((^^.(s)) + (7-7)((/?^(s)))ds. 
Jo Jo 

We verify that the above expression has k — 1 derivatives uniformly bounded since (i) 
X — > is C^~^(X); (ii) 7 has derivatives bounded on X; (iii) {(p^ — <f>x){s) has 

k — 1 derivatives bounded by \\f3 — 

It thus remains to handle the term 

v{x) = / ^{ifxis))ds. 

Jt+{x) 

The function x — > ^{(px{s)) is of class C''~^{X) by regularity of the flow and because 7 is 
of class C^{X). Derivatives of order k — 1 oi v{x) thus involve terms of size i+{x) — t{x) 
and terms of the form 

D^^li^Dt'^'^li^xih)) - t+/^r^^™7fe(t+))), < m < A; - 1. 

Because /3 is of class C^{X), then so is x — 7((^^(s)). Since the latter function has k — 1 
derivatives that are Lipschitz continuous, we thus find that 

This concludes the proof of the proposition when flo = [xq. Applying Lemma 13.81 as 
before, we verify that 

\Dl-' [(/io(x±(x)) - /io(x±(x)))e-/o^'^'^(^^(^»^^] I < C||/io - hWc^^xy 

By the triangle inequality, we deduce the error estimate on /x — /i described in the 
proposition. □ 

Proof of Lemma \3.8[ Let us assume without loss of generality that t^{x) < t^{x). 
We have seen that 

\Mt+{x)) - Vx{t+{x))\ < CU{x)6, 6 = \\P- /3||co(x). 

We also have that \f3{(px{t+)) — l3{ipx{t+))\ < C2S. From Lemma 13.91 we know that 
■ n(x+(x)) > 2Cit+(x) for some constant Ci > 0. 
Let us assume ffist that Cit+ > C2S so that f ■ ei > Cit+(x) where Xq = (fx{t+), 
vq = P{xo), and ei = ?7,(x+(x)). We want to show that the integral curve of /3 starting 
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at Xq at time with velocity Vq exits X in a time of order 5 so that Xj^{x) — is of 

order 5. 

In an appropriate system of coordinates, we have xq = (— ?/, 0) and x+ = (0, c) 
where < y < C6t^. We verify that 



Jo Jo 

'pxo (t) will be outside of the convex domain X as soon as its first component becomes 
non-negative, which implies that 



It remains to choose t+ > AC^C^S to obtain the existence of a time t so that ip^Q {2C36) 
is outside of X. This shows that the distance traveled by ifxo is of size 6 so that 
— x^{x)\ < C6. 

We have treated the case > a6 for some constant a sufficiently large. It 

remains to address the case < a6. For this, a sufficiently positive curvature of dX 

is necessary and we use Lemma 13.91 Indeed, we know that 



We also know that ■ n(x+(x)) > 2aCi6. We can then invoke the second result 

of Lemma [319] and obtain that \x+{x) — x+{x)\ < C6. 

At this stage, we have thus proved that independent of \x+{x) — < CS 

with C a constant independent of 6. The proof of the above result shows that 
\t+{x) — t+{x)\ < C6 as well. Higher order derivatives are now treated in a similar 
fashion. We have seen that 





or equivalently 



Now take t = 20^5 so that the above constraint becomes 



t+ < 4C4C35. 



ip^(t+{x)) - ip:c{t+{x))\ < Ca6^. 



\W{U) - W{U)\ < CU\\D,P - D^PWco^x)- 



Since W{t) is of class and — < C6, we deduce that 



\W{U)-W{U)\<C\\P-P\\c^x) 



which is equivalent to 
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Higher-order derivatives are treated in exactly the same manner providing a bound for 
k — 1 derivatives of — in the uniform norm. 

The error on t+ is obtained as follows. We note that 

After differentiation in space, we obtain 

W{t+{x))Dj+ix) - = - 

Since W is Lipschitz, and t+ — 1+ is small, this implies that 

- = 0(5). 

We have that W is of order and that WifS) = I so that for \i\ sufficiently large, 
W{t^{x)) is invertible. This implies 

iD,U{x)-DMx)) = 0i5). 

Higher-order derivative are treated in the same manner by using the Leibniz product 
rule and the invertibility of W{t+{x)). This concludes the proof of the lemma. □ 

Proof of Lemma \3.9[ Instead of running the characteristics forward from x to 
x^{x), we run the characteristics backwards from x^{x) to X-{x) and show that the 
time spent in X is controlled by the angle the trajectories makes with the normal to X 
at G dX. More precisely, we set y = and v = ipx(t-^{x)) = and 

run characteristics backwards. 

From the equality 

^x{-t) = v - ip^{-s)ds, 







we deduce that 

Wx\-i) -v\< ——, 



and hence 



I / \ / \ I 

\^xK-t) - [y- tv) \ < 



2 



2\t\ 

Let tm the time it takes from to We obviously have that t+(a;) < tm- Let 

By{R) the (unique) ball of radius R tangent to dX at y E dX and such that X C By{R). 

In a system of coordinates with By{R) centered at and v = Ci, we find that 
\y — p = R^ + |f — 2rtv ■ n{y). We want ipx{—t) G X C By{R). This imposes that 

\y-tv\<R + -^^t\ 

Let us define as the maximal time a trajectory spends in X, which is a bounded 
quantity. Then the above imposes that 

\v\H^ - 2Rtv ■ n{y) < ^ ^|m^2 < |^|2^^2^ 
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and in other words that is bounded by ^^^2^_p^ v ■ n(y). 

For the second result, we define Bx{R) as the ball of radius R tangent to dX at x 
and such that X C Bx{R). We again have that \(pxo(f)~{^o+'Vot)\ is bounded by p\vo\H'^. 
In a system of coordinates where n{x) = 62, x = Re2, and Vq = {vq ■ ei)ei + {vq ■ 62)62, 
we obtain that ^Pxoit) ^ X G Bx{R) implies that 



This is equivalent to 



C 

\X + tVo\ < \X - Xo \ + R + TTT^^t^- 



2tRvo ■ 62 + t^\vo\^ < |x - xol^ + ^^t^ + 2\x - xo\R + 2\x - xo\^^ + ^f^- 
This implies that 

(1 - p)\v\H^ < C\x - xol + 2tR\v ■ 62I < C{6'^ + M), 

for some constants C that can be made explicit. Solving this quadratic inequality yields 
that tm is bounded as prescribed. □ 

3.5. ISID with two real-valued measurements 

We are now in a position to state our main stability result in the presence of one well 
chosen complex- valued internal data. We fix 6 and let p = ^ + it"*-. We define u as the 
solution to (I29II with the complex- valued illumination g close to Up\dx- The complex- 
valued internal data are then d = pu. As before, we assume that k>l and yU G C^^^{X) 
so that d G C^^^{X\ C), which implies that (3 and 7 defined in (|53|) are of class C^{X). 
The results of Proposition 13.71 yield the following result: 

Theorem 3.10 Let us assume that {p,q) and {fl,q) are elements in V and that 
\\9 — Up\sx\\co(x) sufficiently small so that u does not vanish on X. Under the 



hypotheses of Proposition 3.1 and assuming that \t\ > |to| with |to| sufficiently large, 
we have that 

Wfi - fl\\ck-i(x) < C\\d - d\\{cHX;C))- (62) 
Moreover, we have the following stability result provided that k > 3: 

Ik - qWc^-Hx) < C\\d - d\\(cHX;C))- (63) 

Proof Let us define /iq = d\Qx/g and fi = d\gx/g. By assumption, g does not 
vanish on dX. We thus deduce that ||/io — Ato||c*(ax) is controlled by \\d — (i||(cfc(j?;C))- 
The rest of the proof of the theorem is a direct consequence of the results obtained in 
the preceding section and of the proof of Theorem 13. 6[ □ 

We now make a few comments on the differences between the 2n— and 2— internal 
data settings. In the presence of n fields, the geometry of X is allowed to be rather 
general. In contrast, the 2— data setting requires much stronger convexity assumptions 
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on X to avoid that integral curves of the vector field be too close to the boundary dX 
for too long, which would result in a severe lack of stability. Since the integral curves of 
the vector fields (5 and /3 are not known a priori, more (Lipschitz) regularity is required 
on the vector fields to ensure that information propagates along near-by trajectories. 
This is the reason for the replacement of k in Theorem 13.61 by "/c — 1" in Theorem 13. 101 



4. Inverse Diffusion with Internal Data (IDID) 

We now return to the diffusion equation with unknown diffusion coefficient D and 
unknown absorption coefficient aa'- 

-V ■ DWu + aau = 0, X, u = g, dX. 

The theorems stated in section [2] are straightforward consequences of the results 
presented in this section in a slightly more general setting. 

Using the standard Liouville change of variables, v = \fDu solves 



At> + gt> = 0, 

with 

A\/D (Ja 



^ D 

The internal data in photoacoustics are given by 

a = = —J=V = jJLV^ jJL 



D VD 

We assume we know V-D on dX. This allows us to prescribe v on dX and thus to 
reconstruct fi and q as in the preceding section. Then we find that 

-A/D-g/D = /i, (64) 



so we can solve for a/D and then get aa = fiVO. 

Let us recall that G F = i/t+^-+2+^(X) C C'+^X). We assume that 

{D,aa) G Ai for some M. This implies that (g, /i) G V for some P. Indeed, that 

solves ( IMl) implies that A + g whose inverse is compact does not have as an 
eigenvalue. Here, as always, k > 1. 

The above calculations show the unique reconstruction of {D, aa) from internal data 
for well-chosen boundary distributions as stated in Theorem 12. 1[ From Theorem 13.61 in 
the 2n— internal data setting, we get the following result. 

Theorem 4.1 Let k > 2 and assume that {D,aa) and {D,aa) are in M. with D\qx = 
D\Qx on dX . Then there is an open set of2n real valued boundary values g in C^'°'{dX) 
for a > \ such that we have the stability estimate 

\\D - DWcki^x) + \Wa - ^aWcHx) < C{M, k)\\d - J||(cfe+i(x))2"- (65) 
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Proof. The main result consists of getting the stabihty on D mentioned above. 
Since k >2, we have stabihty of the reconstruction of q in C''~'^{X) and of /i in C'^(X) 
provided that the boundary conditions are weU-chosen. Thus we have 



(A + q){VD-VD)=i2-jl+{q- q)V D. 



By elhptic regularity, we deduce that {VD — y D) is bounded in C^{X), and hence the 
result. □ 

Finally, we deduce from Theorem 13.101 the following result 
Theorem 4.2 Under the hypotheses of Theorem \3.10i and those of Theorem \4.1 



we 



obtain in the 2-internal data setting the following result. Let k > 3 and {D, 0"^) and 
{D, da) be in M.. Let us assume that D\qx = D\qx- 

Then there is an open set of 2 real-valued boundary conditions g in C^'°'{dX) for 
01 > \ such that we have the stability estimate 

\\D - b\\ck-i(x) + Ika - cra\\c^-i(x) < C{M, k)\\d - d\\ck(X;C)- (66) 
The proof of the above theorem is the same as that of Theorem 14.11 
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